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ABSTRACT 

Published parametric models of the Einstein Cross gravitational lens demonstrate 
that the image geometry can be reproduced by families of models. In particular, the 
slope of the mass-profile for the lens galaxy is unconstrained. However, recent models 
predict a dependence of image flux ratios on the slope of the mass profile. We use 
this dependence to constrain the mass profile by calculating the likelihood of the slope 
using published mid-IR flux ratios (including microlensing variability). We find that 
the galaxy is likely to be flatter than isothermal, and therefore that the mass-to-light 
ratio is decreasing in the inner kpc. 

Key words: quasars - individual: Q22347+0305 - gravitational lensing - microlensing 
- galaxy models 


1 INTRODUCTION 

Q2237+0305 (The Einstein Cross) was discovered in the 
CfA Redshift survey (Huchra et al. 1985). The object com¬ 
prises a quasar at redshift z = 1.695 that is gravitation¬ 
ally lensed by a foreground galaxy (z = 0.0394) produc¬ 
ing 4 images with separations of ~ 1". Many models have 
been proposed for the projected galaxy lens mass distribu¬ 
tion based on observations of the lensed images (e.g. Kent 
& Falco 1988, hereafter KF88; Schneider et al. 1988, here¬ 
after S88; Kochanek 1991, hereafter K91; Rix, Schneider & 
Bachall 1992, hereafter RSB92; Wambsganss & Paczynski 
1994, hereafter WP94; Witt, Mao & Schechter 1995, here¬ 
after WMS95; SWL98; Chae, Turnshek & Khersonsky 1998, 
hereafter CTK98). Studies which employ parametric lens 
models having galaxy and image positions as constraints de¬ 
scribe a degeneracy between the ellipticity of the mass profile 
and its slope v. As a result, the total magnification is virtu¬ 
ally unconstrained. On the other hand, the models of SWL98 
and CTK98 predict variation of the flux ratios with v (this 
dependence arises due to the slightly asymmetric image ge¬ 
ometry). While optical flux ratios cannot be used as model 
constraints due to microlensing and uncertain differential 
extinction (e.g. S88; KF88; K91; WP94), measurements in 
the radio (Falco et al. 1996) should be reliable due to an 
extended emission region. However Q2237+0305 is faint in 
the radio, and the observations have large uncertainties. As 


a result, the majority of published models predict flux ratios 
consistent with those measured in the radio. 

Agol, Jones & Blaes (2000, hereafter AJB00) observed 
Q2237+0305 in the mid-IR and measured flux ratios with 
uncertainties significantly smaller than observations in the 
radio. Each of the 4 images are observed through the bulge of 
the lens galaxy which has an optical depth in stars that is of 
order unity (e.g. KF88; S88; SWL98). This results in a high 
probability for microlensing, and the mid-IR observations 
were coincident with significant microlensing of the optical 
flux. This disparity between optical and mid-IR flux ratios 
has been used to demonstrate that the mid-IR emission re¬ 
gion is large with respect to the microlens Einstein Radius 
(AJB00; Wyithe, Agol & Fluke 2001, hereafter WAF01). 
The large source size indicates that unlike the optical emis¬ 
sion, the mid-IR emission is not subject to significant mi¬ 
crolensing. In addition, extinction is not important in the 
mid-IR. The mid-IR flux ratios should therefore be reliable 
indicators of the true flux-ratios, and can be used to con¬ 
strain galaxy mass models. 

In Sec. [3 we assume microlensing of the mid-IR is neg¬ 
ligible and compute the likeli-hood for different published 
macro-models by comparing the observed mid-IR flux ra¬ 
tios with the predicted flux-ratios. In Sec. we discuss the 
effect of mid-IR microlensing on the results obtained. 
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Table 1. Values for the image C flux at different epochs in the mid-IR (AJBOO, Agol et al. 2001) and V-bands (Wozniak et al. 2000; 
OGLE web page). 


Wavelength 

Date 

Image C flux 

11.7 iim 

28th July & 21st September 1999 

3.15 ±.4 mjy 

V-band 

1st August 1999 

0.46±.01 mjy 

V-band 

26th September 1999 

0.34±,01 mjy 

11.7 fim 

5th October 2000 

3.55±.3 mjy 

V-band 

7th November 2000 

0.180±.015 mjy 


O O 






V 
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Figure 1 . The dependence of magnification //. on power-law index v. The circle represents model 2a from RSB92, the asterix, triangle 
and diamonds represent the models from K91 having internal, mixed and external shear respectively, and squares models from WP94. 
The dotted line joins points for models with zero core radius from CTK98, and the dark lines show values computed from the model of 
SWL98. 


2 CONSTRAINING PUBLISHED MACRO 
MODELS FOR Q2237+0305 

Many authors have attempted to explain the geometry of 
Q2237+0305 via gravitational lensing using various models 
(summarised in SWL88 and WAF01). In this paper we con¬ 
sider the physically motivated parametric model of SWL98 
which includes a bar component to correct for the apparent 
misalignment between the major axis of elliptical power- 
law models and the observed galactic orientation (e.g. K91; 
CTK98). SWL98 found a degeneracy between the ellipticity 


and slope of the mass profile (this degeneracy is a generic 
feature of parametric models (described by WP94)). For 
comparison, we also consider models from four other stud¬ 
ies. Firstly, the model of RSB92 who assumed that mass in 
the lensing galaxy was proportional to observed light. We 
use their model 2a, which did not employ flux ratios as con¬ 
straints. Secondly we consider the model of CTK98 who used 
a tri-axial power-law mass profile with a core (the position 
angle of this model is offset from the galaxy inclination axis 
by ~ 10°). CTK98 incorporated the 3.6cm fluxes measured 
by Falco et al. (1996). However the error-bars on these ratios 
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Figure 2. The dependence of flux ratios R on power-law index v. The circle represents model 2a from RSB92, the asterix, triangle 
and diamonds represent the models from K91 having internal, mixed and external shear respectively, and squares models from WP94. 
The dotted line joins points for models with zero core radius from CTK98, and the dark lines show values computed from the model of 
SWL98. The solid and dashed horizontal lines represent the measurement and 1 <t uncertainty of the mid-IR flux ratios from November 
2000 (Agol et al. 2001). 


were too large to break the degeneracy between ellipticity 
and v. In addition we consider models having a power-law 
mass-distribution with an external shear from WP94, and 
the isothermal models with internal, mixed and external 
shears from K91. 

Fig. |f] shows magnification as a function of projected 
power-law mass-function index v for the models discussed. 
The circle represents the constant mass-to-light model of 
RSB92 (the value of v was inferred from the photomet¬ 
ric deconvolution of Schmidt 96). The asterix, triangle and 
diamonds represent the models from K91 having internal, 
mixed and external shear respectively, and squares models 
from WP94. The dark lines show results computed using 
the model described by SWL98, and the dotted lines join 
points for published results from CTK98 (for a zero core 
radius). The large dependence of the magnification [i on v 
(WP94; SWL98; CTK98) is clearly evident from this figure. 
Six possible flux ratios have been calculated from the 4 im¬ 
ages (Fig. !■ The dependence of R on v is not as strong as 
that of however R is directly observable. 
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Constraints on galaxy models for Q2237+0305 are ob¬ 
tained by comparing the predicted flux ratios with obser¬ 
vations made in the mid-IR. Mid-IR flux ratios are use- 
full for two reasons. Firstly extinction, is not important in 
the mid-IR. Secondly, mid-IR emission is thought to come 
from an extended region, and should not be subject to mi- 
crolensing (Agol, Jones & Blaes 2000). Both these issues 
plauge broad-band optical determinations of the flux-ratios. 
The solid horizontal lines in Fig. |i] are the observed mid- 
IR (November 2000) flux ratios (Agol et al. 2001). Since 
there was no pattern noise in these observations, the errors 
are assumed Gaussian (the horizontal dashed lines in Fig. ^ 
show the lcr errors). There is a large range in predicted flux- 
ratio between different models, however Fig. H shows model 
flux ratios that differ less from the mid-IR observations for 
mass profiles which are shallower than isothermal. Likeli¬ 
hoods for v were constructed for each i mag e ratio, and com¬ 
bined for independent image pairs (Fig. H, left-hand panel). 
The flux ratios predicted by the models of CTK98 suggest 
that an isothermal profile is 3-4 orders of magnitude less 
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Figure 3. Likelihoods for v (left) and ML /,,, i ge / ML sc (right). Two independent likelihoods have been combined in each case. Values 
of likeli-hood smaller than 10 -4 have been plotted at the 10 —4 level. The circle represents model 2a from RSB92, the asterix, triangle 
and diamonds represent the models from K91 having internal, mixed and external shear respectively, and squares models from WP94. 
The dotted line joins points for models with zero core radius from CTK98, and the dark lines show values computed from the model of 
SWL98. 


likely than a mass profile with v = 0.5. Less stringent limits 
are obtained from the models of SWL98. However, a value 
of v = 0.8 is favoured over v = 1.2 by more than a factor 
of 10 2 — 10 3 . Note also that the mid-IR flux ratios push the 
solution towards the deeper part of the x 2 valley shown in 
Fig. 3 of SWL98. The astigmatic models from WP94, the 
constant mass-to-light model from RSB92, and the isother¬ 
mal models of K91 are strongly disfavoured with respect to 
the elliptical power-law models of SWL98 and CTK98. 

While the approach described does not constitute a self 
consistent inclusion of the flux ratios as model constraints, it 
is reasonable for two reasons. Firstly, the possible effects of 
microlensing are difficult to consistently include in a minimi¬ 
sation scheme (see the following section). Furthermore, since 
the relative errors in the 1R fluxes are much larger than the 
errors in position used to constrain published models, con¬ 
sistent inclusion of the Rs as constraints would leave the 
solution in the same x 2 valley in z/-ellipticity space. 

Schmidt (1996) deconvolved HST photometry of the 
lens-galaxy into components of de Vaucouleurs bulge, expo¬ 


nential disc, and Ferrers ellipse bar. At the location of the 
images, most light is contributed by the bulge and disc com¬ 
ponents, and we consider only them (results from Trott 2000 
suggest that the halo component is also negligible at the lo¬ 
cation of the images). Furthermore, the location of the im¬ 
ages coincide with the region where the contributions of the 
bulge and disc components are comparable. Assuming that 
the exponential disc profile can be interpolated into the in¬ 
ner kpc, the slope of the mass distribution can be re-cast as a 
ratio of mass-to-light between the bulge ( ML^uige ) and disc 
( MLdisc ) components. The right hand plot of Fig. [j shows 
the limits re-displayed with v replaced by MLbuige/MLdisc- 
The mass-to-light ratio is likely to be higher in the disc than 
in the bulge. 


3 THE EFFECT OF MICROLENSING 

In this section we consider whether microlensing of the mid- 
IR emission region is an important factor for constraining 
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3.1 


the slope of the lens galaxy from mid-IR flux ratios. Sec. 
discusses the prob ability of the size of the mid-IR emission 
region. In Sec. 3.2 these source-size probabilities, are used to 
compute likeli-hoods for the flux ratios that include possible 
microlensing variability. 


3.1 Limits on the mid-IR source size from two 
epochs of data 

Two epochs of mid-IR observations of Q2237+0305 have 
been made in September 1999 (AJBOO) and November 2000 
(Agol et al. 2001). Both epochs show significant variation in 
flux ratios between the mid-IR observations and the V-band 
observations described by Wozniak et al. (2000). The 1999 
data have previously been used to place limits on the phys¬ 
ical scale of mid-IR emission (AJBOO; WAF01). Briefly, the 
argument is as follows. The optical source is known to un¬ 
dergo significant microlensing (Irwin et al. 1989; Corrigan et 
al. 1991; 0stensen et al. 1996; Wozniak et al. 2000a,b) and 
the variability record has been used to demonstrate that 
the optical source must be significantly smaller than the 
microlens Einstein Radius (e.g. Wambsganss, Paczynski & 
Schneider 1990; Wyithe, Webster, Turner & Mortlock 2000). 
Therefore at any one epoch, the optical flux ratios may dif¬ 
fer significantly form the theoretical mean value (predicted 
from the galaxy mass model). If the mid-IR source is larger 
than the optical region, the mid-IR flux ratios will differ 
from those in the optical, and lie closer to the theoretical 
value. If the mid-IR source is much larger than the microlens 
Einstein Radius, then the mid-IR flux ratios will equal the 
theoretical value. In combination with an assumed macro¬ 
model, the flux-ratios can therefore be used to limit the size 
of the mid-IR emission region. The limits obtained using this 
method are sensitive to the average flux ratio and therefore 
to the macro-model assumed. However a similar argument 
applies to the flux ratios between observations of the same 
image at two different epochs, and this approach is pursued 
below. The flux levels of image C for both the September 
1999 and November 2000 epochs are summarised in Tab. jlj 
The Image C optical flux decreased significantly between 
September 1999 and November 2000, while the mid-IR flux 
remained steady. This temporal variability in the flux-ratio 
can be used to place limits on the mid-IR source size which 
are much less systematically dependent on the macro-model. 
Here we adapt the method described in detail in WAF01 for 
two different image ratios at the same epoch to calculate lim¬ 
its for the same image ratio at two different epochs. Because 
the epochs are separated by only 1 year, a combination of 
transverse velocity and microlens mass (expressed in units 
of the microlens Einstein Radius r/o) must be assumed. Note 
also that differential reddening of the V-band data need not 
be considered since the fluxes are measured from the same 
image. 

We have combined likelihoods {Rir\S ', Sopt) for 
the mid-IR flux ratio Rir given the optical flux ratio Ropt 
(both measured for image C between 1999 and 2000), and 
mid-IR source size Sir (see WAF01) with uniform logarith¬ 
mic Bayesian priors 

Lt'-i prior i 

dlog(SiR) X 

= 0 


where Sir < S max 

otherwise (1) 
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for the unknown mid-IR source size Sir: 

Ps(S < Sir\Sopt) = 

N _C ds ' dJ fF:wr, {R '«' s ’- s °" ) ' ( 2 ) 

N is a normalising constant. Fig. [i] shows the resulting cu¬ 
mulative distributions Ps{Sir < S \ Sopt) for Sir assum¬ 
ing the image C macro parameters calculated by SWL98 
(model xi from WAF01, re = 0.69,7 = 0.71). A uniform 
optical source size was used: Sopt = 0.025 r ] 0 (single pixel). 
The distributions are plotted assuming 4 different upper cut¬ 
offs of the Bayesian prior, Smax = 20,100, 500 and 2500 rj 0 . 
Four cases of epoch separation are shown A r/ 0 = 0.05, 
A rjo = 0.10, A 7] 0 = 0.20 and Ar/ 0 = 2.0. These correspond 
to galactic transverse velocities of ~200, 400, 800 and 8000 

km/sec , where (m) is the average microlens mass. 

The large separation corresponds to an unreasonably large 
velocity or small microlensing mass, and is included for com¬ 
parison as it approximates the case of two random epochs. 
We interpret the 1999 image C peak as a cusp event (Wyithe, 
Webster & Turner 2000). The ~ 1 year duration of this peak 
therefore indicates that a few tenths of an Einstein radius 
separate the 1999 and 2000 observations. Fig. ^ shows lit¬ 
tle variation in the limits obtained for values of separation 
in this range, and we therefore calculate limits at individ¬ 
ual transverse velocities rather than introduce an additional 
prior. 

Because both flux ratios are measured from the same 
image, the average flux ratio equals 1 for all models. The 
limits obtained should therefore be less dependent on macro¬ 
model than those presented in WAF01. Fig. |H] shows cumu¬ 
lative distributions assuming A r/ 0 = 0.1, Sopt = 0.025 7 0 
(single pixel) and the same 4 values of Smax- Distributions 
are shown for each of the macro-models i — xi described 
in WAF01. These models have mean magnifications ranging 
from 2.72-5.56, and the lack of dependence of the distribu¬ 
tions on the macro-model suggests that the mid-IR source 
size limits obtained apply to all reasonable macro-models 
for Q2237+0305. 

We find Sir > 0.5 r/ 0 with greater than 90% ( Sopt = 
0.0257 o ) confidence. As a result of the smaller flux-ratio 
variation, this limit is not as strong as those presented in 
WAF01. Since the mid-IR did not vary by more than the 
observational uncertainty, the data carries no information 
upper limits, which depend on Smax■ Smax is the parameter 
to which the probabilities are most sensitive. 


3.2 model constraints including microlensing 


The results presented in Sec. g assume that the mid-IR flux 
ratios accurately represent the true flux ratios (and hence 
that microl ensi ng of the mid-IR source is negligible). How¬ 
ever in Sec. 3.1, we discussed macro-model insensitive limits 


on the mid-IR source size and it was shown that while this 
should be larger than r/ 0 , there is some probability that the 
source is small enough to be affected by microlensing. In 
this section we consider how microlensing of the mid-IR af¬ 
fects the constraints imposed on the models of SWL98 by 
the mid-IR flux ratios. We have used a ray-tracing code gen¬ 
erously provided by Joachim Wambsganss to compute mag¬ 
nification maps corresponding to the macro-parameters for 
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Figure 4. The cumulative probability for the IR source size P (Sir < S\Sopt) given a uniform optical source: Sopt = 0.02577 o (single 
pixel). The mid-IR source was Gaussian with halfwidth a = Sir/2, and the parameters for image C were k = 0.69, 7 = 0.71. Plots are 
shown assuming 4 cutoffs ( Smax ) on the logarithmic Bayesian prior for Sir ( 20r ) o , solid lines; 100r) o , dashed lines; 500?7o, dot-dashed 
lines; 2500 r) 0 , dotted lines). The 4 panels show limits assuming Aij = 0.05, Arj = 0.10, A 77 = 0.20 and Aij = 2.0. 



Figure 5. The cumulative probability for the IR source size P (Sir < S\Sopt ) given a uniform optical source: Sopt = 0.025ij o (single 
pixel), and a Gaussian mid-IR source with halfwidth a = Sir/2. The macrolensing parameters were for models i — xi described in 
WAF01. Plots are shown assuming 4 cutoffs (Smax) on the logarithmic Bayesian prior for Sir ( 20rj o , solid lines; 100rj o , dashed lines; 
500 770 , dot-dashed lines; 2500 rj 0 , dotted lines). The epoch separation was A 77 = 0.10. 
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Figure 6. Likelihoods for v (left) and MLi u i gf ,/ML,n ac (right) including the contribution of microlensing. Plots are shown assuming 
4 cutoffs ( Smax ) on the logarithmic prior for ,S '/r (20rj o , solid lines; 1 00 //,,. dashed lines; 50(b/ o . dot-dashed lines; 2500?y o , dotted lines). 
Two independent likelihoods have been combined in each case. The macro models were computed according to SWL98. 


each image at 22 different values of v. Three sets of mag¬ 
nification maps were computed to reduce numerical noise. 
Because the distributions of mid-IR source size shown in 
Fig. [ij are not sensitive to the macro-model, they may be 
used to compute a series of source size weighted flux-ratio 
distributions for each v. These distributions have been com¬ 
bined with the mid-IR observations to produce likelihoods 
for v and MLbuige/MLdisc (Figs. |jj). Plots are shown assum¬ 
ing an optical source size Sopt = 0.025 r/ 0 (single pixel), a 
separation of Ar j = 0.1 and the 4 Bayesian priors for Sir 
discussed in Sec. 3.1. The inclusion of microlensing weakens 
the dependence of the likeli-liood on the slope of the mass- 
distribution. It may therefore be said that, in terms of the 
flux ratios the effect of microlensing of the mid-IR emission 
region is important, and this renders self consistent inclusion 
of the mid-IR flux ratios as model constraints difficult. We 
note that the observed mid-IR flux ratios have been used 
both to estimate the source size, and using the probabilities 
obtained, to place limits on v. The two sets of limits are 
therefore not independent. However, the overlap is minimal 
since only the 1999 and 2000 image C mid-IR fluxes were 


used to obtain the source size limits, while all image fluxes 
from 2000 have been used to compute likelihoods for v. 

The left-hand panel in Fig. shows the a-posteriori 
probabilities for v assuming the likelihoods in Fig. |(| and 
a Bayesian prior that is flat between v = 0 and v = 1.3. The 
likelihood appears to flatten below v ~ 0.8 — 1.0, and we 
assume it to be constant below and including v = 0.7. The 
limits from the combination of image ratios A/D and B/C 
rule out a mass density that is steeper than isothermal with 
~ 90% confidence. If the mid-IR emission is due to dust then 
it should be emitted from a region > lpc. We have recom¬ 
puted the results in Figs. [J-[^ using a prior for mid-IR source 
size having a lower cut-off of 10?; o (~ 0.3 , < ' rn ' 1 pc). The exclu- 

Lm q 

sion of the possibility for small sources limits the microlens¬ 
ing variability, and results in stronger limits of u < 0.8 with 
90% confidence. 


4 CONCLUSION 

Mid-IR flux ratios (Agol et al. 2001) have been used to con¬ 
strain degeneracies in the slope the mass profile for published 
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Figure 7. A-posteriori probability for v (including the contribution of microlensing). A flat Bayesian prior is assumed between v = 0 and 
v = 1.3. The likelihoods are assumed constant below and including v = 0.7. Plots are shown assuming 4 cutoffs ( Smax ) on the logarithmic 
prior for Sir (20 r] 0 , solid lines; 100r; o . dashed lines; 500?7 o , dot-dashed lines; 2500 r/ 0 , dotted lines). Two independent likelihoods have 
been combined in each case. The macro model was described in SWL98. Two figures are shown. The left-hand nlots are computed from 
the likeli-hoods in Fig. [(j. The right hand plots show the corresponding limits where the calculations in Figs, [ij-y have been re-computed 
using a lower cut-off on the prior for mid-IR source size of 10// o . 


models (Chae, Turnshek & Khersonsky 1998; Kochanek 
1991; Rix, Schneider & Bahcall 1992; Schmidt, Webster & 
Lewis 1998; Wambsganss & Paczynski 1994) of the bulge of 
the lensing galaxy G2237+0305. Microlensing of the mid- 
IR is at a much lower level than for optical emission but 
may also be important. For the model of Schmidt, Webster 
& Lewis 1998 our calculation includes microlensing of the 
mid-IR emission by first finding the a-posteriori probabil¬ 
ity for source size in a macro-model insensitive way using 
the variation of the image C V-band to mid-IR flux ratio 
between 2000 (Agol et al. 2001) and 1999 (Agol, Jones & 
Blaes 2000). Microlensing induced flux ratio distributions 
were produced by combining distributions for the mid-IR 
source size with magnification distributions corresponding 
to macro-models having different power law slopes. Likeli¬ 
hoods for the slope of the mass profile were then computed 
from the combination of mid-IR flux ratio observations (with 
Gaussian errors) and microlensing induced flux ratio distri¬ 
butions. The mid-IR flux ratios favour a mass profile which 


is flatter than isothermal, resulting in an cllipticity lower 
than the light profile. Following the photometric disc-bulge 
deconvolution of Schmidt (1996), this can be interpreted as 
a mass-to-light ratio which is higher in the galactic disc than 
in the bulge. A shallow mass profile results in high magni¬ 
fication, and this helps alleviate the difficulty of reconciling 
the observed quasar luminosity with microlensing imposed 
source size constraints (Rauch & Blandford 1991). 

Future mid-IR monitoring with reduced uncertain¬ 
ties will improve constraints on lens galaxy models for 
2237+0305 in two ways. Firstly, if no variability of the mid- 
IR is observed, limits on the mid-IR source size will become 
tighter, and microlensing a less important consideration. On 
the other hand, if variability is observed there will no longer 
be strong dependence of mid-IR source size limits on the up¬ 
per cutoff of Bayesian prior. Secondly, at the current level 
of uncertainty, the observations differ by only ~ 3<r from 
the predicted flux ratios for the least favoured models (from 
Schmidt, Webster & Lewis 1998). Reduced uncertainties will 
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therefore increase the disparity in the likelihoods of different 
models and result in tighter constraints 
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